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ABSTRACT 


Transient performance of gas turbines has a strong bearing on output and 
component life. For this reason, several articles have been written on the 
dynamic simulation of gas turbine systems in electrical generation, 
cogeneration, and marine applications. These models provide a basis for 
this present work. This paper describes a mathematical and computer model 
that was developed to investigate the dynamic response of a simple (no 
reheat, regeneration, or other auxiliary equipment) single-shaft gas turbine 
system. The model uses design parameters normally incorporated in gas 
turbine design (e.g. load coefficient, flow coefficient, and deHaller Number) 
as well as compressor and turbine stage geometry and compressor and 
turbine material properties. Also incorporated is a combustion chamber 
model. Other input parameters are included to enable the model to be 


adaptable to various system sizes and environments. 


The results of several trials are displayed, showing how given parameters 
(e.g. fuel-air ratio, load, and efficiencies) were varied and results were 
compared. A thermodynamic balance resulted in a 26% system thermal 
efficiency for the reference state parameters. Also, a start-up procedure 


similar to an industrial gas turbine system was simulated with this model. 





Examination of the results from the trial runs show the computer model to be 
reliable, given that no control system was incorporated. Because of the 
approach used, the computer model is easily adaptable to further 


improvements and a control system. 
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Hub-to-Tip Ratio 

Compressor Blade Height (m) 
Compressor Disk Radius (m) 
specific Heat (J/kg K) 
Fuel-to-Air Ratio 

deHaller Number 

Enthalpy (J/kg) 

Moment of Inertia (kg m*) 


Ratio of Specific Heats 


Mass Flow Rate (kg/s) 


Pressure (Pa) 


Power (kW) 

Radius 
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Temperature (K) 
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Efficiency 
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Torque (N m) 

Flow Coefficient 
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Subscripts 
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Axial 

Compressor 

Load 

Combustion mixture (fuel and air) 
Root 

stage 

Tip 

Turbine 
Stagnation 
Compressor Inlet 
Compressor Outlet 
Turbine Inlet 


Turbine Outlet 
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INFRODUCTION 


Transient performance of gas turbines has a strong bearing on output and 
component life. Issues of thermal fatigue [1] can influence the start-up 
sequence, since rapid start-ups can impair turbine blade life. As hot gas is 
generated in the combustion chamber from a cold-start, gradual (as opposed 
to sudden) heating of the turbine components reduces thermal stresses 
associated with the start-up. Consequently, the turbine is activated by 
external power with air inlet guide vanes (IGVs) closed, while the combustion 
chamber is simply in a heating duty for about ten minutes. Synchronization 
with the electrical load in the case of power generation also requires 


accurate speed control under varying load. 


In marine propulsion applications, system response in terms of fuel 
consumption for an oscillatory load is also of concern [2]. Different control 
strategies were analyzed in this work, namely Proportional Integral and 
Linear Quadratic Regulator. Both were shown to offer strengths and 
drawbacks with regards to operational aspects. This work shows the great 
value that dynamic modeling can offer in terms of system evaluation prior to 
actual construction. Inevitably, the conclusions depend on the model 
adopted for the gas turbine, as described in [3]. The model used here for 


the gas turbine has a definitive empirical slant, since it is largely based on 
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performance data and time lags acquired from an actual gas turbine test bed. 


This approach was expanded and refined in [4]. 


Rather complex systems have been modeled based on similar empirical 
approaches [5]. In this case, a cogeneration plant consisting of a gas 
turbine, a boiler, and two steam turbines (high and low pressure), were 
simulated. The system was approximated by using the performance curves 
given by manufacturers, and adding necessary time lags based on 
manufacturer information. Time-dependent equations were integrated with 
Euler’s method. It was concluded that dynamic simulation was a valuable tool 


for assessing the suitability of controls and safety measures. 


A more basic approach to modeling was adopted in [6]. The basic 
thermodynamic equations for mass and energy balances were used for 
model formulation, as opposed to manufacturer's curves. Yet, the inertia of 
rotors and shafts was absent from the equations, which probably explains the 
need found by the authors to add time lags to this model. The model 
effectively predicted behavior during load perturbations at full load. A slightly 
more recent work focused on developing the gas turbine model based on 
design parameters [7]. Load coefficients and stage efficiencies were 
employed for this simulation, along with part-load losses extracted from the 


literature. Again, inertial aspects seem to have been neglected by the author. 





The need to combine basic thermodynamic equations with inertial aspects for 
generic models was recognized recently. In [8], the transient equations for 
inlet duct pressure and mass flow rate were formulated. Compressor and 
combustion chamber conditions were calculated from a transient analysis 
coupled with polytropic process equations. The thermal inertia of all 
components was neglected, and performance maps of each component were 
used to obtain isentropic efficiencies. A dynamic equilibrium equation was 
used to calculate the acceleration of the turbine. Good comparisons with 


experimental performance were obtained. 


The present work is inserted in this context: whereas most models are 
semiempirical, in that they rely on time constants and performance maps, 
serious recent attempts to extend the generality of gas turbine models have 
been made. Clearly, the flexibility of a model to competently handle a wide 
variety of turbines is the target. Similar to previous work, the trend of 
formulating the basic mechanical and thermal equations for each component 
is followed here. In an extension of previous work, it is attempted wherever 
possible to include mechanical and thermal design parameters to broaden 
the reach of the model. Here, the main contribution consists of modeling the 
compressor stage by stage, and of including thermal and mechanical inertial 


effects based on first-principles, rather than empirical constants. 





It is realized that complex machines cannot always be represented by simple 
models. Hence, the limitations of this approach are clear, especially at part- 
load. However, with the advent of computing power and of transient 
simulation software, it is inefficient to restrict modeling to numerical fitting of 
performance curves. Based on a literature search and this present work, it 
is ascertained that models of this type will be used more and more in the 
future to establish feasibility of design and of control and operating 


strategies. 





DEVELOPMENT OF THE MATH MODEL 


The mathematical model of the turbine system is based upon the three main 
components as seen in Figure 1 below: compressor, combustion chamber, 
and turbine. In addition to these main components, there is a component of 


mass flow rate included in the model. 


COMBUSTION 
FUEL CHAMBER 





COMPRESSOR TURBINE 


Figure 1. Schematic of the three main components: 
compressor, combustion chamber, and turbine. 


Assumptions 


These assumptions are used in the development of the gas turbine model: 


(1) For the purposes of inertial calculations, compressor and 
turbine stages can be modeled as disks with rectangular 
blades attached on the edge of the disk. 
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(2) 


(3) 
(4) 
() 
(6) 


(7) 
(8) 


(9) 


(10) 


The increase in size of each successive stage of the 
compressor and the decrease in size of each successive 
stage of the turbine are approximated by constant factors. 
Rotors in the compressor and turbine are perfectly balanced. 
The compression and expansion processes are adiabatic. 


Air acts as an ideal gas. 


The compressor inlet pressure is equal to the turbine outlet 
pressure. 


There are 16 compressor stages and 4 turbine stages. 


Stage efficiencies for the compressor and turbine are 
constant. 


The combustors can be modeled as open-ended cylinders of 
constant wall thickness. 


The temperature drop across the turbine is divided equally 
among the four stages. This distribution may be altered as 
desired. 


The turbine and compressor, along with the load, are connected to a single 


shaft, such that the dynamic equation of the shaft is: 


Each component of the model is described below. 
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Compressor and Turbine Stages 


To model the moments of inertia of the rotating masses in the compressor 
and the turbine, each stage is modeled as a disk with rectangular blades 


attached as seen in Figure 2a: 





Figure 2a. Model of compressor and turbine stage. 


This moment of inertia model is derived from the model of a rotating disk 
found in reference [9]. Although the cross-sectional free area is a function of 
the design flow rate and blockage coefficients, a simpler approach was used 
to model the radial dimensions. First, a cross-section of the stage radius was 


assumed, as seen in Figure 2b. The dimensions shown are relative to the 
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design radius of the stage (r). Each of the cross-hatched areas was modeled 
as a rotating ring, and the section at the center line was modeled as a 
rotating disk. These sections were added for the total moment of inertia for 
one stage, and then modeled as a single rotating disk, as previously seen in 
Figure 2a. The blades were then added into the moment of inertia 


calculation. 
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Figure 2b. Cross-section of a compressor and turbine stage. 
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For each successive stage in the compressor (decreasing in size) and each 
successive stage in the turbine (increasing in size), a constant “decrease 


factor’ (compressor) and “increase factor’ (turbine) was used to estimate the 





change in size (see Figures 3a and 3b below). The mass of all 16 
compressor stages were summed for a total mass, and the same was done 
for the 4 turbine stages. Using these parameters, the total moment of inertia 


was calculated for the compressor (I-) and the turbine (I+). 





Figure 3a. The compressor is modeled Figure 3b. The turbine is modeled 
as a series of 16 approximately linearly as a series of 4 approximately linearly 
decreasing stage sizes. increasing stage sizes. 





Compressor 


The model of this component is based upon the specification of design 
variables normally employed to characterize a given turbine design. Hence, 
the compressor model was developed using load coefficient, deHaller 
Number, flow coefficient, degree of reaction, and stage geometry. Starting 


with degree of reaction and load coefficient, the equations used are [10]: 


y = hee (2) 
Ghpeat,) 

Fy ee 3 
No» ES No, 


and for the blade tip speed, U, 
U = (radius + blade_ height) - Q (4) 


Combining equations (2), (3) and (4), the temperature rise in each 


compressor stage can be calculated as follows: 


(cdr + cbh)° - y - Q? 
Toss = Tors + 7 io Oa (5) 


Ncs ° Mp 


where “cdr” is the compressor disk radius and “cbh’ is the compressor blade 


height. 
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successively using equation (5) for each stage results in the compressor exit 
temperature. The pressure rise in each stage is found by using the following 


equation: 





(6) 


The compressor isentropic efficiency, nc, is found using the relationship 


between polytropic and isentropic efficiencies [11]: 


(k-1) 


Ss 
P, 


Noe (kK) (7) 





The torque required by the compressor is found by again using degree of 
reaction, deHaller Number, load coefficient, flow coefficient, and stage blade 


tip speed: 


h, —h, 


aT (8) 
oz — Ao 


where (Nop - Nox) Is: 


Ve V2 
Noo — Ag, = h, a Me) = hh se *) (9) 





Combining equations (8) and (9) and letting 


Noo = No = Ah, 


results In 


Substituting in the deHaller number, H, 


H = wa (10) 
my 
results in 
1 Z 2 
R 1 Se) 11 
Definitions of Flow coefficient and Load coefficient are: 
V 
= — 12 
bag (12) 


Ah 
Y = ie (1 3) 








Taking equations (12) and (13), and substituting them into equation (11) 
results in a relationship between degree of reaction, load coefficient, and flow 


coefficient: 


1 2 2 
pete eet} 
R = 14+ ———___—_ (14) 
WV 


Ah, = 2-U-(V,, - V,.) = We (15) 


and combining this with the load coefficient defined as [10]: 


results in 
W 
= le (17) 


Using two definitions of power, 
P=Q-1, and P=mw, 
equating and solving for the stage work, Ws, results in: 


W, = Q + Ts 





m 





Combining this with equation (17), 


yeas (18) 





Equating equations (16) and (18), 
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The stage torque can be found as 


medal = eye ay 


ce 


Additionally, the compressor efficiency must be taken into account. Recalling 


the load coefficient, equation (13), 





(ideal) 


=a (actual) 





Using this efficiency, and the fact that the temperature differences are 
Summed across all the stages, the total compressor torque can be calculated 


as follows: 


to = ee an 08) (20) 
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Turbine 


The torque imparted from the fluid onto the turbine is developed starting with 


the work of an individual stage of the turbine: 
We = Comic - (Ts - Ta)g (21) 


Substituting this into the load coefficient definition [10]: 


2-W, 
y= (22) 
and from equation (4), 
U = (tdr + tbh) - © (23) 


where “tdr” is the turbine disk radius and “tbh” is the turbine blade height. 
Substituting equation (23) into equation (22), 


— 2 We (24) 
Y  ~dr + tbh)’ - Q? 


Again using the definitions of power, 





and solving for ts: 


W, -m 
Q 


Coa = 





Using the fact that the mass flow rate through the turbine is the mass flow 


rate of air plus the mass flow rate of fuel, so the stage torque becomes 


Si — (25) 


The turbine efficiency is factored into the above expression as follows: 


Ah 
y= ie (ideal) 


_ Ah, 
 mgleieen 





(actual) 


Where the turbine efficiency is found using the relationship between 


polytropic and isentropic efficiencies [1 1]: 


(k-1)-nts 
k 


1 (28) 





In addition, the turbine outlet temperature is calculated and knowing the 


ambient pressure and turbine inlet temperature as follows: 





Sy 
T, = Ty (EB (27) 


Substituting equation (21) into equation (25) and incorporating the efficiency 


results in the turbine torque equation: 


ria eee Cee (ToT, 
ee (28) 


The load coefficient of the turbine is calculated using the temperature drop 


across the turbine as a parameter: 





(29) 


C, (ATs } 


= 30 
‘ 2 - (tdr + tbh)’ - Q? - n, (89) 





Mass Flow Rate 


An equation for mass flow rate through the compressor is given by [11]: 


Z 
® r 
apap lt-(E) ly, (31) 
t 


Making the following substitutions into equation (31) 


r = cdr + cbh 


F, 
r 


Vv, = 6-2 


the mass flow rate can be calculated as follows: 


m = p-n- (cdr + cbh)’ -(1- B2)- 6-2 (32) 


where the density is calculated using the ideal gas law. 
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Combustion Chamber 


The model for the combustion chamber is developed using the thermal mass 
of the chamber itself, the thermal mass of the fuel-air mixture, and the 
thermal mass of the first stage of the turbine. A control volume analysis of 


the combustion chamber reveals the following differential equation: 


(me) [Cp - Ty | + (m, | . (LHV) 
z [ms |C, . Te + (mC, ) - Tos (33) 
where mCp is the thermal mass of the system, which is modeled as follows: 


mC, = [p- Cp - Voll +[p- Cp - Volfuc-ar + [p-Cp - Volfis (34) 
stage 


The pressure drop across the combustion chamber is calculated by: 
P, = P, + (1-8) -(P, -P,) (35) 


where “8B” is the fractional pressure drop across the combustion chamber. 





pi 


The complete mathematical model is comprised of equation (1) as the basic 
equation, using equation (20) for tc, equation (28) for tz, calculated values 
for I- and ly, and user-input load and friction torques. The independent 
variable in these equations is time (t). The remaining equations in the 


mathematical model are: 


massmiow tate iano iecceeonamts. cs (32) 
compressor outlet temperature ................. cee eeeeeee eee (5) 
compressor outlet PreSSUre........... eee ec eeceeceeeee ees (6) 
compressor isentropic efficiency .................c ccc eececeeeee ee (7) 
turbine isentropic EFFICIENCY ............. ccc ecccseececesseeeeenee (26) 
turbine outlet temperature.................. ec eceeeeeeece eee es (27) 
turbine load coefficient.....................ccccccssescesscecceureeeses (30) 
turbine inlet temperature..................... (33) and (34) 


combustion chamber pressure drop.................000..000ee (35) 
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DEVELOPMENT OF THE COMPUTER MODEL 


Similar work has been done in simulating gas turbine systems. References 
[3] and [4] describe a linear model, while [8] describes a non-linear model 
similar to this paper. However, this model provides inputs for compressor 
and turbine stage geometry and design parameters such as B (hub-to-tip 
ratio), y (load coefficient), and o (flow coefficient). Additionally, this program 
integrates a combustion chamber model. Several other inputs are required 
to run the model, including ambient conditions, compressor and turbine 
efficiencies, and compressor and turbine material properties. These input 
parameters enable this model to be adapted to a wide variety of gas turbine 
systems and environments. It is important to note that this model does not 
incorporate a control system as in [2] and [12]. However, the model is easily 


adaptable to such a system. 


The software package VisSim ® was chosen for this model because it 
provides “a visual environment for developing continuous, discrete, multirate, 
and hybrid system models and performing dynamic simulations [13].” There 
are several other special features about VisSim ® outlined in [13] that make 
this a superb platform. The program itself is called “gtsim1.vsm” created by 


the author. 
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Figure 4 is the “top level” of gtsim1.vsm which shows the basic logic and 
interconnections of the previously discussed mathematical model. The 
program incorporates equation (1) as the basic equation, using equation (20) 
for tc, equation (28) for t7, calculated values for Ic and I7, and input load and 
friction torques. The load and friction torques are computed as constant 


factors multiplied by the square of the angular velocity [sec '}’. 


Additionally, gtsim1.vsm uses the following equations for the associated 


parameters: 
ila SS 1lOW FAIS .. ::c., 2. see oe 5: eich sSa eens «see (32) 
compressor outlet temperature ................ cece eeec eee eee (5) 
compressor outlet PreSSure............. ccc cecceeceeseee eens (6) 
compressor isentropic effIClENCy.......... ee cecccceeece seca ee (7) 
turbine isentropic efficiency ............... cece ceeceeeeeeee ees (26) 
turbine outlet temperature.............. eee eee eee eee eens (27) 
turbine load coefficient... cee eecceeceeeeeeeeeues (30) 
turbine inlet temperature .............0000.0c (33) and (34) 


combustion chamber pressure drop......................:00000 (35) 





24 


‘pessedoid Buyjaq ,,s{euB)s,, JO U0}}99Ijp BY) QJed|PU] SpEBYMOY 
‘wesBoid ey Jo 2j60) 91seq oy Buymoys ‘wWsa’}uUjs76 Jo JeAe] GOL ‘py eunBi4 




















SNOSINWdWNOO a ~SiINdNI|] _— -. % 
GNV S3AYNND SLAdNI OINWNAG Se | ll SLNNI 
JONVWHOSY3d Y¥3HLO -OWNSHL| | YOLSNAaWOD) [| SNIGYNL || YOSSSYdWOD 
s nS Agen seinbue 7 a ae oe —— Ayooja seinbue™ 
; J * -CRYG . BIVOUI! JO as 
[Avie Sermon Yl | a. ~ | juawou euIquN) 
sejnBue | . 
7 aap lal 
NU! JO 
JUBWOW JOSSAIGWOO 
—— [eo fees . 
ri es 
"peo von | SONOMA 
“4S i | peo! }« le - jas . 
a ee . ae 
. z Meek he rs 3Nbj0O} peo} 
= ou 
= 1 
~— anbso) eulqun) (11; vba lq 
1 ..——)- ——--—___-,.-—-. ob. [g 
= ame poe Ce = See 3 enbio) JOSsasdWw0o | | 
ae, ae Ol | Me cig ee 0 ee i 
P - kyooyan seinBue SSWW Pee woyu i | 
GOU819}4IP 
eunyesodwea)} 
ROR yanNopgiul 
ae | ee zs x a mre s - eS ; euiqiny 
ae mp ee sue | ee oe ee ee | sNolWinoWwo 
Zh Zi JNIGYHNL 
el 3 s 
> 


P kyon seinBue Ayoojan seinBue 





25 


RESULTS 


To demonstrate the versatility of this model, several parameters were varied 
independently. In each case, all variables were held to their reference state 


(Appendix A) except the variable under discussion. 


Reference State Run 


Figure 5 shows the adopted baseline run. This run is based on a start-up 
from “cold” conditions, at constant fuel/air ratio, with a frictional load 
proportional to the square of the rotational velocity. Note that in this plot, the 
angular velocity is left in units of [s"'] in order to clearly fit into the plot scale. 
The model allows an understanding of the variables in the uncontrolled start- 
up scenario. As long as the turbine produces no work, T3 (peak cycle 
temperature) and T, (discharge temperature) are quite close. The 
compressor discharge pressure increases gradually with time. The air- 
breathing nature of this turbomachine is clearly reflected here: as soon as 
the speed increases and the air mass flow ramps up, compressor pressure 
and peak temperature start to increase sharply. As both T3 and P2 increase 


toward the design value, the turbine power output also increases. 
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Congruently, T, decreases, and the mass flow rate rapidly increases towards 


its design value. At steady-state conditions, the turbine develops 56.6 MW to 


drive the load. 
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Figure 5. Response of turbine to baseline conditions (fuel/air ratio = 0.013, standard 
atmospheric conditions, and friction load and shaft load proportional to angular 
velocity squared). Start-up is from “cold” conditions at constant fuel/air ratio with 
frictional load. Inlet guide vanes (IGVs) are completely open. 
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Energy Balance 


A control-volume energy balance was performed on each component of the 
model to ensure realistic parameters were passed between components and 
realistic output power values were obtained. Table 1 below gives values for 
this analysis. For this system at the adopted baseline parameters (Appendix 


A), the overall thermal efficiency was found to be 26%. 





Component Power Out Component Power In 
(MW) (MW) 


es 
ee | 
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Table 1. Thermodynamic balance of the adopted baseline run. 
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Start-Up Runs 


In the following figures and tables, a given parameter was varied and the 
system allowed to start up from initial conditions. The resulting steady-state 


speed and time to reach the steady-state speed for each case is shown. 


In Figure 6, the fuel-air ratio (held constant from turbine start-up to steady- 
state) was varied with the same start-up rotational velocity of 30 RPM. Up to 
approximately 500 seconds, the speed is the same for all three cases (f/a = 
0.012, fla = 0.013, and fla = 0.014). After approximately 500 seconds, it is 
seen that for greater fuel-to-air ratios, less time is needed to reach higher 


Steady-state speeds. 


—t— fla = 0.012 ' 


_—O— fla = 0.013 | 


| —X— fla = 0.014 | 


RPM 
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Figure 6. Influence of fuel to air ratio on response time and final velocity. 
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Figure 7 demonstrates the effect of compressor polytropic (small stage) 
efficiency [11] on steady-state speed. The plot shows that the rotational 
speed of the system increases as the compressor polytropic efficiency 
increases. This is due to the ability of a more efficient compressor to impart 
more work to a given amount of air. The greater work results in increased 


peak pressure, which in turn increases turbine output. 
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Figure 7. Influence of compressor polytropic efficiency on response time and 
rotational velocity . 
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Figure 8 demonstrates the effect of turbine polytropic (small stage) efficiency 
[11] on steady-state speed. The plot shows that as the turbine polytropic 
efficiency increases, the rotational speed of the system increases. This Is 
due to the ability of a more efficient turbine to take more work from the 
combustion gases and use it to drive the compressor, friction load, and shaft 


load. 
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Figure 8. influence of turbine polytropic efficiency on response time and 
rotational velocity. 
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Tables 2 and 3 below demonstrate the results of varying compressor load 
coefficient ( y ) and compressor blade height, respectively. Here in Table 2, 
increasing the load coefficient means that more work is communicated per 
stage to the air. Hence, the peak pressure increases, resulting in a more 


efficient cycle which leads to greater power output. 


Load Time (sec) to 
Coefficient Steady-State Steady-State 
y RPM RPM 
0.25 3296 925 
0.30 3450 725 
Uso 3536 615 


Table 2. The effect of varying compressor load coefficient ( y ) (the 
ratio of stage work to blade kinetic energy). As the load coefficient is 
increased, the final speed also increases. This is due to the 
increased amount of work per stage as speed increases. 
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In Table 3 below, it is seen that the model responds reasonably as the blade 
height is increased. The longer blades can sweep more air at a given 


rotational velocity, resulting in increased power output. 


Compressor Time (sec) to 
Blade Height Steady-State Steady-State 
(m) RPM RPM 
0.39 3256 950 
0.41 3450 Zo 
0.43 3629 aro 


Table 3. The effect of varying compressor blade height. As the 
blade height is increased, there is more “swept” area through which 
more mass can flow during a revolution of the compressor stage. 
The direct result is that for more mass flow, more power, and hence, 
greater steady-state speed. 
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Variations During Steady-State Operation 


For these runs, the model was allowed to ramp-up to its reference steady- 
state conditions from a “cold” start. The amount of time taken to get to these 
reference conditions was approximately 725 seconds. From the steady-state 
conditions, certain parameters were varied as described below and results 


compared. 


Figure 9 shows the effect of load changes during steady-state operation of 
the gas turbine. “Load - Decreased” indicates that the constant load factor 
(the factor multiplied by the square of the angular velocity to obtain the actual 
load) was instantaneously reduced from 1.2 to 1.0. “Load - constant’ 
indicates no change in the load, and “Load - Increased” indicates that the 
constant factor was instantaneously increased from 1.2 to 1.4. The load was 
decreased and increased at different times for clarity. The results are as 
expected: an increased load results in lower rotational speed and a 


decreased load results in higher rotational speed. 
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Figure 9. Influence of changes in load during steady-state operations. 
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Figure 10 is similar to Figure 9, except that the fuel-air ratio instead of the 
load is changed during steady-state operations. In Figure 10, “fla 
decreased” indicates that the fuel-air ratio was decreased from 0.013 to 
0.012, “fla constant’ indicates no change, and “f/a increased indicates that 
the fuel-air ratio was increased from 0.013 to 0.014. Again, the parameter 
(fuel-air ratio) was decreased and increased at slightly different times for 
clarity. The results are again as expected: increased fuel-air ratio results in 
higher rotational speed and a decreased fuel-air ratio results in lower 


rotational speed. 
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Figure 10. Influence of changes in fuel-air ratio during steady-state operations. 
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Figure 11 shows the effect of closing the air inlet guide vanes (IGVs) to the 
compressor. At steady-state conditions, the IGVs are 100% open, allowing 
the system to ramp-up to its reference mass flow rate (338.7 kg/s). In Figure 
11, the “constant” line shows the IGVs remaining at 100% open. The other 
two lines represent, as noted, |GVs closed to 50% and completely closed 
(0%). The plot shows that the effect of closing the IGVs is reduced mass flow 


through the compressor, hence, less power output of the system. 
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Figure 11. influence of changes in mass flow (closing inlet guide vanes) during 
steady-state operations 
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Start-Up Routine 


Start-up procedures of an actual industrial gas turbine system are given in 
[12]. These procedures are facilitated using an automatic control system and 
an external “start-up” device such as a diesel engine or an electric motor. 
Despite the lack of a control system, an investigation was done on 


gtsim1.vsm to see if these start-up procedures could be simulated. 


To facilitate this test, a “start-up” engine was simulated by changing the sign 
of the load term in equation (1), making it negative. The mass flow rate and 
fuel input were set to zero. Using a VisSim ® “slider” block, the load was 
manually varied from a more negative value to a less negative value, 
allowing the system to come to 100% of design (reference state) rotational 


speed. Figure 12 was generated from this run: 
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Figure 12. Trial run of start-up procedures given in reference [12]. 


The plot of T, does not match reference [12] due to the lack of fuel input 
during this run. In reference [12], there was fuel input, giving the system 
outlet temperate higher values than in Figure 12. This start-up trial of 
gtsimi.vsm did not exactly simulate reference [12]; however, it does 
demonstrate that this model can be adapted to an external start-up device 


and control system to simulate that in reference [12] or others. 
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CONCLUSIONS 


A mathematical and computer model was developed to investigate the 
dynamic response of a simple (no reheat, regeneration, or other auxiliary 
equipment) single-shaft gas turbine system. The model used parameters 
normally incorporated in gas turbine design as well as compressor and 
turbine stage geometries and compressor and turbine material properties. 
Also incorporated was a combustion chamber model. Input parameters were 
included to enable the model to be adaptable to various system sizes and 


environments. 


The results of several trials were displayed, showing how given parameters 
(e.g. fuel-air ratio, load, and efficiencies) were varied and results were 
compared. A thermodynamic balance of the adopted baseline state resulted 
in a 26% system thermal efficiency. Further, a literature search showed that 
industrial gas turbines are brought on-line from “cold” conditions via an 
elaborate start-up procedure. The goal of this procedure is dual: to avoid 
undue thermal stresses leading to fatigue, and to allow ready synchronization 
to external loads. The turbine model presented in this paper can respond to 
similar requirements as imposed manually by the programmer. In an 


uncontrolled start-up, the turbine does converge quite rapidly to steady-state 
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conditions. It is noteworthy that the peak temperature increases gradually at 
first, and more rapidly after the air flow is established. Most of the 


temperature increase is quite gradual, on the order of 50K/min. 


Examination of the results from the trial runs show the computer model, 
gtsim1.vsm, as a functional gas turbine model, given the fact that no control 
system was incorporated. For a frictional load and constant fuel/air ratio, the 
turbine model converges to a steady-state condition. In this condition, the 
turbine and internal and external loads are in balance. Small departures of 
load and fuel/air ratio from steady-state values result in small departures in 
rotational speed. Because of the approach used, the computer model is 


easily adaptable to further improvements and a control system. 
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APPENDIX A 


Reference State Parameters 


Compressor Variables 
Material Density sense eee exec ceeresrcscxerese-.-0areeeasease-s 7000 kg/m” 
BIG Supe IG Neer cote aco. oacn0cs8e2seccacsedenchestlowts.csls ce 0.41 m 
BB Slee UE eae cs iiss. candenkds.ccccssecceeseeelaeee 0.02 m 
BISGGeMIMIGKMCSS .wrrmmentet st s).2o) 00. 04.cce.widscvsleevecusesseseeseeeeee 0.004 m 
PIS I AC US occa ck sce cantscaossoeccastcsecenesisucoseeereeneeee 0.52 m 
[DISK IAICKMNESS scien eiitact sein ae sie ivdaeas odedeetaddoasenas sae 0.02 m 
Nit MDE ROM DIOO SS sc me we eat cat Sc. cites: we dsdacs neseteeee so eesa eee anes 62 
Decrease Factor (of successive stage SIZ@S) ...........00 ee, 0.98 
POlVILODIG ill ClOINey meetin 8.52... cot secs) vc. Jodsdanesdiae atlas ee eee 0.85 


Turbine Variables 


METGITE) BYE cess acaissne a eee eee. 7000 kg/m? 
Material Specific Heat......................ccccscsccnssesccesecenecsssanuenes 582 J/kgK 
PAG St eC UG ence eee es 5 okocs Uitnsicesevecesaa cee ee 0.5m 
Blade ube Ctl ates) en i ae inact asks 0.02 m 
BIAS [LNIGKNESS, smn trnis ss ccciccvscsesdeete  eeN no annccosansceeeeenes 0.004 m 
PES Gig e ACMMES <r I ov Rleiggalsa'k cbs os'dGccain Ee 0.4m 
DISKO MICKMESS - 22ers OMY coo. scaan's sch gonemegesneeeeneeee 0.02 m 
NLIMNID CG Gtits) AC © SINT oes oss. hw suadais saean ne ee eee onc. ee 8/ 
Decrease Factor (of successive stage SIZ@€S) .............c eee ee a 
FOU NG tenn [te CY cotter oars dca sc ccs ogee esas. eee ee Sse 0.89 


Combustor Variables 


Combustor Diameter ..0........... 0. cece ccc ececececcececcccecececseececeeeceees 0.35 m 
Gomisustor Lemaitre: see. ..........0.002c.0c0see ste ences vevesnucecness: 0.7m 
Number of CombuStofs 2000.00.00... ccc c ccc ccccccececceccceccceccccceccceccccesceceeces 8 
Combustor Wall Thickness ........0..... 0. cece cce cece cc eccccecceccceceseeees 0.008 m 
Combustor Material Density .........0...0..cccccccccccecsscessseeseeeee 4510 kg/m” 
Combustor Material Specific Heat..........000.0cc eee. 582 J/kgK 


FRESSULS POI pee es a 3 Sec yak ive phee esd cncshehee Hi eaeazeee 2% 
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Thermodynamic Variables 


[UMBIRESSULC. .. ccna teenet reset tcctees «- ccs sacaieescssccsvansieauaeeoente 101325 Pa 
POU ICs SPO MMOLE Os. Nee ac ea saa odes ccs e. bees ce nndadanccesoaaoSeeeees. 289 K 
RatiorOh Specific MEAS TOMA -.........05....0.....ccssssvreececes sone MMlMscteess. 1.4 
Ratio of Specific Heats for combustion mixture ..................0...0.. 1.33 
Universal Gas Constant...................ccccccccecceeseceeseeeeeeeens 8314 J/molK 
SPCCHiGumeae Of All... cckec: sess 2ess.ce..cesadvoumametnacsasdaas 1005 J/kgkK 
Specific Heat of Combustion Mixture ...................ceeee 1100 J/kgK 
Lower Heating Value of Fuel...........0..... eee eens 90000000 J/kg 
Molecular Weight of Alr...................ccccccssecesesceeesereeees 28.85 kg/kmole 
Other Variables 
External Load Factorl..........:..::ccccscceccceccecccecccesccacceceescceeeeeees 2.0 kgm 
Load Coefficient (Compressor) ..............ccc ccc cceccsceecceeeccueseueseuunees 0.3 
Flow Coefficient (Compressor and Turbine) .......................ceee ee 0.72 
Hub-to-Tip Ratio (COMPIeSSOLS)............... cece ceecceeeece secu eeeeeeeeeeees 0.75 
Frictiompload) F ACtOr MMe << occ cvececese ccs esseccveseuttereess- este 0.2 kgm? 


Results of gtsim1.vsm using these reference parameters 


Fuel-Aik Ratio (INPUut).. GeRe....0..5....cccccccssesceecesseasnoncaeecsesees linc. 0.013 
Initial RPM Condition on Integrator Block ......................6e. 30 RPM 
Turbine Inlet Temperature ................ccc cee ccc ceceeeceeceeeeeeeseeees 1265 K 
Turbine Outlet Temperature.............0.cc ccc ee ceneeeeeeeeeneeuees 680 K 
Rotational VelOCILY ..... Mets ss vecscvevccverowsvsescocae MMe rxes teases 3450 RPM 
ERESSUGC Rt One seer PP AN Sir nas os csastwesavenauacmemeimetiiecs-ia.0essacs 17.4 
Time to reach steady State............. eee eee approx. 7/25 seconds 
WaSS RIOW IRAte <.. coca ooo sauces oes sessccbemendessesseesees 338.7 kg/s 
POWeR OUI OUT si ccccccres MIMI reeea acces. aan pete sssnssccacdgaaMaeesssssessaaueh 56.6 MW 
Compressor Isentropic Efficiency..........0..... ccc cccccceeeesseeeseueeees 0.78 


Turbine Isentropic Efficiency... ccc cecceecceeceeseeesuseeeeees 0.93 





Appendix B - Survey of Turbomachinery Handbook/1995 [14] 








Shaft Turbine 

Speed Pressure Inlet Exhaust Exhaust 
Company Model (RPM) Ratio Temp (K) Temp (K) Flow (kg/s) 
Mechanical Drive Turbines 
ABB GT-35 3450 12.0 1123 649 92.0 
European RLM2500 3600 18.4 1501 797 69.0 
European RLM2500+ 3600 21.0 1501 785 83.0 
European RLM5000 3600 25.4 1443 107 121.0 
Kvaerner LM2500-PE 3600 18.4 1086 801 68.0 
Kvaerner LM6000-PA 3600 29.6 1126 130 125.7 
Electrical Generation Turbines 
ABB GT26 3000 30.0 1508 893 543.0 
ABB GT24 3600 30.0 1508 883 376.0 
ABB GT13E2 3000 14.6 1373 797 532.0 
ABB GT13D 3000 11.9 1263 763 394.0 
ABB GT11N2 3600 15.0 1358 797 375.0 
ABB GT11N 3600 Vo 1300 778 317.0 
Ansaldo V84.2 3600 10.8 sa 821 353.0 
Ansaldo V84.3 3600 16.1 ea 823 433.0 
Ansaldo V94.2 3000 11.0 eos 816 514.0 
Ansaldo V94.3 3000 oO. ait 823 624.0 
Dresser DR-63G 3600 29.6 1104 736 172.0 
Mitsubishi MW501 3600 14.0 ane 794 305.0 
Mitsubishi MW701 3000 14.0 sia 786 453.0 
Mitsubishi 501F 3600 14.6 aa 853 437.0 
Mitsubishi 701F 3000 16.0 bei 826 667.0 
Mitsubishi 501G 3600 19.0 a 866 544.0 
MTU LM2500+ 3600 23.0 ee 799 83.0 
MTU LM2500+ 3000 23.0 a4 784 84.0 
Siemens V84.2 3600 10 i 817 360.0 
Siemens V84.3A 3600 16.6 =n 835 454.0 
Siemens V94.2 3000 11.1 ia 813 519.0 
Siemens V94.3A 3000 16.6 ee 835 540.0 
Marine Drive Turbines 
FiatAvio LM2500 3600 19.2 1130 826 69.8 
GEC LM2500 3600 18.0 aa 813 70.0 
Kvaerner LM2500 3600 18.4 1086 801 68.0 
MAN FT8 3600 19.8 . 1131 84.1 
MTU LM2500+ 3600 23.0 on 799 83.0 
Pratt& Wh FT8 3600 18.7 une 716 78.1 
Rolls Royce WR21 3600 127 a 629 713.0 
Turbo Power FT8 3600 19.8 ce 716 83.9 
Average Values .............. 3429 18.1 1294 801 286.5 

Shaft Pressure Turbine Exhaust Exhaust 


Speed Ratio Inlet 
(RPM) Temp (K) 


Temp (K) Flow (kg/s) 
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